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Abstract. The radiative return offers the unique possibility for a measurement of the cross section of 
electron-positron annihilation into hadrons over a wide range of energies. The large luminosity of present 
4>- and B-factories easily compensates for the additional factor of a due to the emission of a hard photon. 
Final states with photons at large angles can be easily identified. The rate for events with collinear photons, 
however, is enhanced by a large logarithm and allows, in particular at lower energies, for a complementary 
measurement. The Monte Carlo generator PHOKHARA, which includes next to leading order corrections 
from virtual and real photon emission, has been extended from large photon angles into the collinear re- 
gion, using recent results for the virtual corrections. In addition, the present version includes final state 
radiation for muon and pion pair production and final states with four pions. Implications for the experi- 
mental analysis at three typical energies, 1.02, 4 and 10.6 GeV, are presented: the magnitude of these new 
corrections is studied, possibilities for the separation of initial and final state radiation are proposed, and 
the differences with respect to the previous treatment based on structure functions are investigated. 
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1 Introduction 

Measurements of the cross section for electron-positron 
annihilation into hadrons are essential for the interpreta- 
tion of the recent, precise results for the muon anomalous 
magnetic moment a M Jl]. Similarly they are relevant for 
our knowledge of the running of the fine structure constant 
and thus crucial for the analysis of electroweak precision 
measurements at high energy colliders [3EE]- 

Of particular importance for these two applications is 
the low energy region, say from threshold up to centre-of- 
mass system (cms) energies of approximately 3 GeV for 
a M and 10 GeV for a(Mz)- Recent measurements based 
on energy scans between 2 and 5 GeV have improved the 



* Work supported in part by BMBF under grant number 
05HT9VKB0, EC 5-th Framework EURIDICE network project 
HPRN-CT2002-00311, EC 5-th Framework contract HPRN- 
CT-2000-00149, TARI project HPRI-CT-1999-00088, Polish 
State Committee for Scientific Research (KBN) under contract 
No. 2 P03B 017 24, MCyT Plan Nacional I+D+I (Spain) un- 
der grant BFM2002-00568 and Generalitat Valenciana under 
grant CTIDIB/2002/24. 

a e-mail: czyz@us . edu . pi 

b Supported in part by 'Marie Curie Training Site' at Karl- 
sruhe Univ.; e-mail: grzel@joy.phys.us.edu.pl 

c e-mail: jk@particle.uni-karlsruhe.de 

d Supported in part by E.U. TMR grant HPMF-CT-2000- 
00989; e-mail: german.rodrigo@cern.ch 



accuracy in part of this range [H]- Similar, or even further 
improvements below 2 GeV would be highly welcome. The 
region between 1.4 GeV and 2 GeV, in particular, is poorly 
studied and no collider will cover this region in the near 
future. Improvements on the precise measurements of the 
pion form factor in the low energy region by the CMD2 
and DM2 collaborations [S] , or even an independent cross- 
check, would be extremely useful, in particular in view 
of the disagreement between e + e~data and the analysis 
based on isospin-breaking-corrected r decays 0] . 

Traditionally the energy dependence of the cross sec- 
tion was deduced from experiments, where the beam en- 
ergy was varied over the range dictated by the energy 
reach of the collider. This 'energy scan' allows, at a first 
glance, a fairly simple interpretation of the measurement 
in terms of the so called R-ratio, which enters the afore- 
mentioned applications. Nevertheless also in this case ini- 
tial state radiative corrections (ISR) give rise to complica- 
tions and require a complicated unfolding procedure dis- 
cussed below. 

As an alternative the 'radiative return' has been sug- 
gested [7H51|S1| 10| as a particularly attractive opt ion for o- 
and B-meson factories. These collider experiments oper- 
ate at fixed energies, albeit with enormous luminosities. 
BABAR and BELLE at 10.6 GeV, CLEO-C in the region 
between 3 and 5 GeV and KLOE at 1.02 GeV are the 
experiments of interest for the subsequent considerations. 
This peculiar feature of a 'factory' allows the use of the 
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radiative return, i.e. the reaction 

e + (p 1 ) + e-(p 2 )^ 7 (fci)+7*(Q)Hhadrons) , (1) 

to explore a wide range of Q 2 in a single experiment. 

Nominally invariant masses of the hadronic system be- 
tween the production threshold of the respective channel 
and the cms energy of the experiment are accessible. In 
practice, to clearly identify the reaction, it is useful to con- 
sider only events with a hard photon-tagged or untagged- 
which lowers the energy significantly. 

To arrive at reliable predictions for differential and for 
partially integrated cross sections, including kinematical 
cuts as used in experiments, a Monte Carlo generator is 
indispensable. The inclusion of radiative corrections in the 
generator and in the experimental analysis is essential for 
the precise extraction of the cross section. For hadronic 
states with invariant masses below 2 to 3 GeV, it is desir- 
able to simulate the individual channels with two, three 
and up to six mesons, i.e. pions, kaons, ry's, etc., which 
requires a fairly detailed parametrisation of various form 
factors. 

A first program called EVA was constructed some time 
ago |2] to simulate the production of a pion pair together 
with a hard photon. It includes initial state radiation, final 
state radiation (FSR), their interference, and the domi- 
nant radiative corrections from additional collincar radi- 
ation through structure function (SF) technique A 
similar program that simulates the production of four pi- 
ons together with a hard photon has been developed in 
0. More recently a new Monte Carlo generator called 
PHOKHARA was developed. It includes, in contrast 
to the former generators, the complete next-to-leading or- 
der (NLO) radiative corrections. 

The first version of PHOKHARA incorporates ISR 
only and is limited to ir + Tr~j('j) and /i + /i~7(7) as final 
states. PHOKHARA exhibits, however, a modular struc- 
ture that simplifies the implementation of additional ha- 
dronic modes or the replacement of the current (s) of the 
existing modes. Its first version was designed to simulate 
configurations with photons emitted at relatively large an- 
gles, 9 2 » m 2 e /s. In this case it is legitimate to drop terms 
proportional to m 2 , an assumption that leads to a consid- 
erable simplification of the virtual corrections ^2] ■ Subse- 
quently analytical results for the virtual corrections, that 
are also valid into the small angle region, were obtained 
in The extension of the program PHOKHARA into 
this small angle region, incorporating these new analytic 
results are the central topic of the present paper. The de- 
scription of this new feature and numerous tests of the 
program stability and technical precision are contained in 
Section 2. 

Final state radiation can affect the measurement of the 
pion form factor, and quite generally of the R-ratio. How- 
ever, using suitable cuts, its effects can be significantly 
reduced. Moreover, given sufficiently large event rates its 
magnitude can be extracted experimentally by varying the 
cuts and/or comparing events with different photon angles 
with respect to beam and pion directions, respectively. 
The charge asymmetry that arises from ISR-FSR inter- 



ference provides another important handle on this 'back- 
ground'. For this reason FSR from 7r + 7r~ and has 
been incorporated in the upgrade of the program and will 
be discussed in Section 3. The final state is still 

limited to its QED part, e.g. the narrow resonances (J/V') 
are not (yet) included. 

The program has also been extended to include final 
states with four pions, following the lines discussed in 9 . 
The implementation of these new channels will be dis- 
cussed in Section 4. 

2 The radiative return for small-angle 
emission and tests of the program 

The study of events with photons emitted under both large 
and small angles, and thus at a significantly enhanced 
rate, is particularly attractive for the 7r + 7r _ final state 
with its clear signature |15Ul(Slll7lf l8.19.20 . In contrast 
events with a tagged photon, emitted at a large angle, 
have a clear signature particularly suited to the analysis 
of hadronic states of higher multiplicities |21U22| . 

The inclusion of radiative corrections is essential for 
the precise extraction of the cross section, which is nec- 
essarily based on a Monte Carlo simulation. The com- 
plete NLO corrections have recently been implemented in 
the program PHOKHARA. However, just like the earlier 
EVA, this program was designed for photon emission at 
large angles ('tagged photons'). For nearly collinear pho- 
tons, corrections from virtual and real photon emission, as 
well as Born terms, must include those contributions pro- 
portional to m 2 /s and even to m^/s 2 , which are enhanced 
by their highly singular angular dependence and thus in- 
tegrate to terms of order 1 and proportional to a/n for 
Born and NLO terms respectively. These mass-suppressed 
terms are significantly smaller than the leading, logarith- 
mically enhanced pieces; they must nevertheless be taken 
into account for a consistent treatment. The evaluation 
of corrections from virtual plus soft photon emission to 
reaction valid for the full angular region, has been 
treated in [2]. Essentially it consists of the calculation 
of the leptonic tensor i^, which has to be multiplied by 
the hadronic tensor H^, so that a fully differential dis- 
tribution is obtained. To arrive at a reasonably compact, 
numerically stable result, the limit m 2 /s -C 1 for is 
considered. However, terms proportional to m 2 must be 
kept if these exhibit the singular angular dependence dis- 
cussed above. 

The differential rate for the virtual and soft QED cor- 
rections is thus cast into the product of a leptonic and 
a hadronic tensor and the corresponding factorised phase 
space: 

da = ■^-L liv dl$>- 1 {$>\ +P2;Q,h) 

dO 2 

x H» v d$ n {Q; qi ,--- ,q n )-%- , (2) 

where d<P n (Q; qi, ■ ■ ■ ,q n ) denotes the hadronic n-body 
phase space, including all statistical factors, Q 2 is the in- 
variant mass of the hadronic system and d<p2(pi J rP2', Q,k\) 
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Fig. 1. Comparison of the Q 2 differential distribution for different values of the soft photon cutoff: w = 1CP 3 vs. 1CP 4 and 
w — 10~ 4 vs. 10~ 5 , at t/s = 1.02 GeV. One of the photons was required to have energy > 10 MeV. No further cuts were 
applied. 
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Fig. 2. Comparison of the Q 2 differential distribution for different values of the soft photon cutoff: w = 10 -3 vs. 10 -4 and 
w = 10 -4 vs. 10 _J , at ^fs = 10.6 GeV. One of the photons was required to have energy > 100 MeV. No further cuts were 
applied. 



is the two-body phase space of the photon and the hadronic 
system. The tensor depends on the four- vectors p±, 
P2, Q, ki and the soft photon cutoff w = E^^/^fs. Its 
explicit functional form is given in The description 
of the hadronic system is model-dependent. It enters only 
through the hadronic tensor 

H »v = jv r \ i (3) 

where the hadronic current has to be parametrised through 
form factors 023,24 25 . The running of a is not taken 
into account in this program and can be included trivially 
in the final experimental analysis. 

The matrix element for the emission from the initial 
state of two real hard photons, i.e. E lt > Wy/s, with i = 



Table 1. Total cross section (nb) for the process e + e _ — > 
7r + 7r~7 at NLO for different values of the soft photon cutoff. 
Only initial state radiation. One of the photons with energy 
larger than 10 MeV for ^/s = 1.02 GeV and larger than 100 
MeV for yjl = 10.6 GeV. Q 2 < lGeV. No further cuts applied. 

w 7s =1.02 GeV 10.6 GeV 

10" a 36.999 (3) 0.15557(7) 

10 -4 37.021 (3) 0.15548(6) 

10~ 5 37.021 (3) 0.15545(6) 



1,2, 

e+{p 1 ) + e-{ P2 )^ 1 *{Q)+ 1 {k 1 )+ 1 {k 2 ) , (4) 
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is calculated numerically following the helicity amplitude 
method with the conventions introduced in |26U27| . The 
results from |12| . which were used for tagged photon events 
are equally applicable for the present purpose. 
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Fig. 3. Comparison of the virtual+soft and hard contributions 
to the 7r + 7r - differential cross section with inclusive analytical 
results. Soft photon cutoff: w = 10 -4 . 



The virtual plus soft photon contribution and the hard 
one depend separately on the soft photon cutoff w used to 
regulate the infrared divergences of the virtual diagrams. 
The former shows a logarithmic ^-dependence. The latter, 
after numerical integration of the phase space, exhibits the 
same behaviour, whereas their sum must be independent 
of w. To show that this indeed occurs is therefore a ba- 
sic test of the performance of the program. The value of 
w that optimises the efficiency of the event generation, 
avoiding at the same time the appearance of the negative 
weights, is determined by this procedure. 

Table ^ presents the total cross section for radiative 
production of a pair of pions calculated for several values 
of the soft photon cutoff at two different cms energies. The 
energy of one of the photons was required to be larger than 
10 MeV for ,/s = 1.02 GeV and larger than 100 MeV for 
a/s = 10.6 GeV. No further kinematical cuts were applied, 
thus allowing to test in particular the small photon angle 
region. 

For i/i = 1.02 GeV the comparison between w = 10~ 3 
and 10 -4 indicates a residual w-dependence. The excellent 
agreement between w = 10~ 4 and 10 -5 , within the error 
of the numerical integration, confirms the w-independence 
of the result. A value around w = 10 -4 seems to be the 



best choice as observed before for large angle photons |12l 

EBl 

In Figs.2] an d0the Q 2 dependence of the differential 
cross section da/dQ 2 is compared for different choices of 
the cutoff, after integration over the remaining kinematic 
variables. Again for = 1.02 GeV the comparison be- 
tween w = 10~ 3 and 10~ 4 shows a residual w dependence 
(Fig. 2^), which disappears beyond w = 10~ 4 (Fig. QJd). 
At a cms energy of 10.6 GeV the result is numerically sta- 
ble for w = 10~ 3 already (Fig. Et)- Stable results are also 
obtained for w around and below 10~ 4 (Fig. Thus 
w — 10~ 4 is used as the default value in the program. 
Similar tests were performed for the four-pion channels 
(see Section 4). 

The present implementation of PHOKHARA covers 
the full angular region for photon emission. This allows for 
a number of tests and comparisons with analytical results 
that were not possible with the previous version. In Fig. 
the results of the program are compared with the analyti- 
cal results from Ref. We use their Eqs. (2.25)+(2.26) 
for the virtual plus soft part and Eq. (2.28) for the hard 
emission part. As it was necessary to change several cou- 
plings in the original formulae, we repeat (Appendix B) 
the expressions actually used for the comparison. Agree- 
ment within the statistical uncertainty and in any case 
better than 10~ 3 is evident from this comparison. 

Initial state radiation is dominated by photons at small 
angles. Inclusion of events with nearly collincar photons 
thus leads to a significant enhancement of the observed 
event rate. The comparison between the differential cross 
sections for large angle photon events (30° < 6 1 < 150°) 
and without restriction on Q 1 is shown in Fig. The pion 
angles are always assumed to be restricted to the region 
30° < 9 7r ± < 150°. Results are presented for two different 
cms energies ( y/s — 1.02 and 10.6 GeV). There is a big 
quantitative difference between these two energies. For y/s 
= 1.02 GeV a huge contribution from small angle photons 
is observed for the full range of Q 2 . In contrast the gain in 
the cross section for y/s — 10.6 GeV is small as a conse- 
quence of the conflicting kinematical constraints of small 
photon and large pion angles. FSR has not been included 
in these figures. For 10.6 GeV its contribution is negligi- 
ble, while for 1.02 GeV and for the cuts used for Fig. 0] 
it is sizeable, but can be reduced using the cuts discussed 
below. 

The R-ratio can in principle be deduced either from 
the measurement of the hadronic cross section, which re- 
quires a precise control of the luminosity, or from the ratio 
between hadronic and event rates. Various radiative 

corrections, e.g. from the running of the electromagnetic 
coupling and from ISR cancel in the ratio between the 
hadronic and rates. Indeed, one obtains by con- 

struction unity, if one considers the properly normalised 
ratio 



4(1 + 2m 2 l /Q 2 )p„ 

ft I K I 2 



■7(7) 



(5) 



where 
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Fig. 4. Comparison of the Q 2 differential cross sections for sfs = 1.02 (left) and 10.6 GeV (right). The pion angles are restricted 
to 30° < d w ± < 150°, while the photon(s) angles are not restricted (upper curves) and restricted to 30° < 9 1 < 150° (lower 
curves) . 




and is the pion form factor. 

The result p^^ = 1 is independent of the restrictions 
on the photon angular region and is true in Born and NLO 
approximations. The phase space of hadronic and 
final states, however, must be fully integrated (Fig. |5Ji). 
For realistic cuts on pion and muon angles the ratio devi- 
ates significantly from 1, a consequence of their markedly 
different angular distributions. The size of this effect de- 
pends on the details of the cuts on photon, pion and muon 
angles as demonstrated in Figs. 03 andj^t- In both figures, 
one observes a significant, few per cent, difference between 
Born and NLO predictions for p n)Ji , depending on the de- 
tails of the cuts on the photon and charged particle angles. 
At 10.6 GeV the ratio p^ is of course again equal to 1 if 
pions and muons are fully integrated (Fig.UJi). In contrast 
to the situation at lower energies, the inclusion of realistic 
cuts does not alter this picture drastically, a consequence 
of the high correlation between photon and pion or muon 
angles: photon and charged particles are essentially emit- 
ted back to back (Fig. Et>) ■ 

The shape of these curves depends only on the pion 
and muon angular distributions, but not on the form fac- 
tor itself. These results can thus be directly used to de- 
duce efficiencies of specific experimental cuts in a model- 
independent way, since pion and muon angular distribu- 
tions are fixed by general considerations. For more com- 
plicated final states (e.g. 4-7T, KKn, ■ ■ •) the correspond- 
ing ratio would, instead of |F^| 2 /3^/4, directly involve the 
corresponding R-ratio, if no cuts on the hadrons are ap- 
plied. Otherwise the results depend on the model for the 
hadronic form factor implemented in the program. An im- 
portant advantage of the radiative return is implicit in all 



these considerations: by measuring Q 2 directly, the invari- 
ant squared mass of the hadronic final state, one has direct 
access to R at the corresponding value of Q 2 . This differs 
from the measurement of the inclusive cross section as a 
function of y/s (energy scan). To extract the true R(s), an 
unfolding has to be performed, which requires in princi- 
ple the knowledge of the cross section over the full energy 
range below and a precise knowledge of the radiator func- 
tion. In contrast, using the radiative return method, it is 
still necessary to know the QED radiator function, but 
no unfolding is required as one measures the Q 2 of the 
hadronic system, and thus has 'access' to the hadronic 
cross section at that given Q 2 . 

Let us discuss those am 2 ./ s and a(m 2 /s) 2 terms, which 
in combination with their singular angular dependence in- 
tegrate to corrections of order a/ir. These are included in 
the present version of PHOKHARA. The corresponding 
leading order corrections proportional to m 2 are typically 
of the order of a few per cent |31|. while the non-leading 
ones are of order 0.1%. This can be seen from Fig. The 
size of these effects is consistent with the expectations for 
a 1% terms without logarithmic enhancement. These terms 
will become important when the precision of the measure- 
ment will be below 1%. Their proper treatment is in that 
case crucial as they do depend on Q 2 and do change the 
Q 2 distribution from which the hadronic cross section is 
extracted. 



3 Initial versus final state radiation 

A potential complication for the measurement of the pion 
form factor or generally of the R-ratio may arise from the 
interplay between photons from ISR and FSR. Their rel- 
ative strength is strongly dependent on the photon angle 
relative to the beam and the pion directions, the cms en- 
ergy of the reaction and the invariant mass of the hadronic 
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system. FSR from hadronic final states cannot be pre- 
dicted from first principles and thus has to be modelled. 
The model amplitude can nevertheless be tested by consid- 
ering charge-asymmetric differential distributions, which 
arise from the interference between ISR and FSR ampli- 
tudes 7 . In leading order the complete matrix element 
squared is given by 



|M| 2 = |Misr| 2 + |Mfsr| 2 + 2Re[Mi SR M^ SR ] 



(6) 



which is still independent of the model for FSR. 

FSR and its interference with ISR were already in- 
cluded in EVA [7] for the two-pion case. The pions were 
assumed to be point-like, and scalar QED was applied to 
simulate photon emission off the charged pions. ft was 
demonstrated there that ISR dominates for suitably cho- 
sen final states, namely those with hard photons at small 
angles relative to the beam, well separated from the pi- 
ons. FSR can therefore be reduced to a reasonable limit, 
and moreover, can be controlled by the simulation (see 
also |32|h Similar results can be obtained using the new 
version of PHOKHARA, were FSR and ISR-FSR interfer- 
ence are included for two pions and muons (see Appendix 
A for details) at LO. The photon emission from pions is 
again modelled by a point-like pion-photon interaction. 
The upper two curves of Fig. |S] describe the differential 
cross section with arbitrary photon and pion angles. The 
contribution of FSR is clearly visible. Once photon emis- 
sion is restricted to angles close to the beam and if the 
pion- and photon-allowed angular ranges do not overlap 
(lower curve), FSR is clearly negligible. 

The third term in the right-hand side of Eq. JHJ, ISR- 
FSR interference, is odd under charge conjugation, and 
its contribution vanishes after angular integration. It gives 
rise, however, to a relatively large charge asymmetry and, 
correspondingly, to a forward-backward asymmetry 



A(9) = 



N T ' (9) 



N* + (9) + N* + (n 



(7) 



The asymmetry can be used for calibration of the FSR am- 
plitude, and fits to the angular distribution A{9) can test 
details of its model dependence. Given sufficiently large 
event rates this procedure can be performed for different 
9-y, thus allowing for an unambiguous reconstruction of 
the FSR amplitude. 

This is illustrated in Figs. and where the an- 
gular distributions of ir + and li + respectively are shown 
for different kinematical cuts. The angles are defined with 
respect to the incoming positron. If no angular cut is 
applied, the angular distribution in both cases is highly 
asymmetric as a consequence of the ISR-FSR interfer- 
ence contribution. If cuts suitable to suppress FSR, and 
therefore the ISR-FSR interference, are applied, the dis- 
tributions become symmetric. 

These investigations can also be performed for different 
photon energies, thus exploring FSR in different regions of 
Q 2 . We will return to this aspect in a future publication. 

At B-factories, where one has to deal with very hard 
tagged photons, the kinematic separation between the pho- 
ton and the hadrons becomes very clear. For events where 
hadrons and photon are produced mainly back to back, 
the suppression of FSR is naturally accomplished and no 
special angular cuts are therefore needed to control FSR 
versus ISR at higher energies (Fig.lll[). The relative size of 
the FSR is of the order of a few per mil, but does depend 
on the value of the pion form factor at y/s = 10 GeV, 
which is extrapolated from the low energy data. 

The suppression of FSR contributions to 7r + 7r~7 events 
is also a consequence of the rapid decrease of the form 
factor above ~ 1 GeV. It is therefore instructive to study 
the corresponding distributions for /i + /i~7 final states. For 
a/Q 2 < 1 GeV FSR is still tiny. Around 3 GeV a small 
charge-asymmetric interference term becomes visible (Fig. 
1120 . which is still irrelevant after averaging over /i + and 
fi~. At large Q 2 , however, FSR plays an important role 
both for the charge-asymmetry and the charge symmetric 
term. 
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4 The four-pion mode 

Because of the modular structure of PHOKHARA, ad- 
ditional hadronic modes can be easily implemented. The 
four-pion channels (27r + 27r~ and 27r°7r + 7r~), which give 
the dominant contribution to the hadronic cross section 
in the region from 1 to 2 GeV, are a new feature of our 
event generator. 

Isospin invariance relates the amplitudes of the e + e~ — ► 
27T + 27r - and e + e~ — ► 27r°7r + 7r~ processes and those for t 
decays into ir~3ir° and n + 2iT~ir 33 9. The description 
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Fig. 13. Diagrams contributing to the 4n hadronic current. 



of the four-pion hadronic current follows [§ll24| . The basic 
building blocks of this current are schematically depicted 
in Fig. El and described in detail in [5j . 

Results obtained with PHOKHARA for these channels 
have been compared with the Monte Carlo, which simu- 
lates the same process at LO ||] and includes additional 
collinear radiation through the SF technique. Typically, 
differences of order 1% are found (see Figs. 1141 and 115(1 . 
which are of the expected size and of the same order as 
for the two-pion final state \12\. 

The generation of the pion four momenta is however 
different from the one described in 9| . In the present ver- 
sion of the program we absorb the most prominent peaks 
in the four-pion hadronic current to obtain a more efficient 
Monte Carlo generation. The Q 2 distribution is peaked 
around y^Q 2 = 1.5 GeV, with a large width of ~ 0.5 
GeV. This is the result of an interplay between several 
resonances present in that region. Nevertheless one Breit- 
Wigner resonance provides an adequate approximation for 
efficient generation. For the approximant and the genera- 
tion of the Q 2 distribution we use 



MQ 2 ) 



- Q 2 (Q 2 - m 2 ) 2 + r 2 m 2 
with m = 1.5 GeV and r = 0.5 GeV. 



(8) 
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It takes care of soft photon emission (s — Q 2 ~ E 1 ) 
and the aforementioned resonant behaviour. For the pro- 
cess e + e _ — > 27r°7r + 7r-7(7) we furthermore absorb the 

to peaks in the four momentum squared Qqh = = 

(Po + P+ + P-) 2 - T ne approximant used for that pur- 
pose, according to which the three-particle four momenta 
squared are generated, reads 



Table 2. Total cross section (nb) for the process e + e~ — > 4-7T7 
at NLO for different values of the soft photon cutoff at y/s = 
1.02 GeV. Only initial state radiation. One of the photons with 
energy > 10 MeV. Q 2 < lGeV. No further cuts applied. 

w 2n + 2n~ 2n°TT + n' 

~ UT 3 0.170167(15) 0.55725(5) 
10" 4 0.170413(14) 0.55844(5) 



10" 



0.170431(15) 0.55845(5) 



i ;4; Q234 



)=2 



r 



(Q 



2 

134 



,2\2 



r 2 ,r 



(Q 



2 

234 



,2\2 



) 2 + r 2 ml 



(9) 



where Q134 and Q234 are the four momenta of the two 
ir°TT + Tr~ subsystems. The other variables are generated as 
described in 

Also for the four-pion modes, we have tested the in- 
dependence of the result of the soft-hard separation pa- 
rameter w. The results are very much similar to the it + tt~ 
mode and are summarised in Tableland Figs. llt)l and ll7l 
Again the choice w — 10~ 4 is the proper one. 



5 Conclusions 



The Monte Carlo generator PHOKHARA, which simu- 
lates the radiative return at electron-positron colliders, 
has been extended from large angles into the collinear 
region using recent results for the virtual corrections to 
photon emission, which are valid for all photon angles. 
Comparing the program with analytical results, a tech- 
nical precision better than 0.5 x 10~ 3 is demonstrated. 
The importance of NLO corrections for the extraction of 
a correct value for the R-ratio is emphasised. 

A number of corrections vanish if the ratio between 
hadron and muon pair cross section is considered. For low 
energies, around 1 GeV, the ratio depends strongly on the 
cuts on the charged particles and corrections have to be 
applied. At higher energies, around 10 GeV, and for low 
Q 2 , the dependence on these cuts is drastically reduced. 

In the new version of PHOKHARA, described in this 
article, final state radiation in leading order treatment is 
included. We discuss the implications for the measurement 
of the pion form factor. Suitable cuts allow, on the one 
hand, the determination of this model-dependent ampli- 
tude, on the other hand it is possible to select configura- 
tions that are entirely dominated by initial state radiation. 

Finally we extend the program final states with four- 
pions configuration as a first step towards the inclusion of 
a multitude of exclusive states. 
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A The implementation of the final state 
emission 

Final state emission of one photon and the final-initial 
state interference terms are implemented in the present 



program in lowest order by means of the helicity ampli- 
tude method for both 7r + 7r~7 and final states. The 
notation is the same as in an d will not be repeated. 
The pion-photon interaction is adopted from scalar elec- 
trodynamics. 

The helicity amplitudes describing the initial emission 
read 



Misr(A 6 + , A e - , Ai) = —qT~ I 

uJ(pi,A e +)-A u I (p 2 ,\ e -) + v\ I (p 1 ,X e +)B unfa, A e -)| , 

(10) 
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where 



and 



.4 



and 



B 



(e*(fci,Ai)-fc+-2e*(fc 1 ,Ai)-pi) J~ 
2fa - pi 

J- (2e*{fa,X 1 )-p 2 -k+e*{fa,X 1 )-) 
2ki ■ p 2 



(e*(fci,Ai) + fcr-2e*(fci,Ai)-pi) J+ 
2ki -pi 

J+ (2e*(fa,X 1 )-p 2 -k^e*(fa,X 1 )+) 



11) 



2ki ■ p 2 

The current J^ 1 for 7r + 7r~ in the final state reads 

J^=ieF 2v (Q 2 ) (q« + -q„-r , 
while for fi + fi~ in the final state it is given by 



(12) 



(13) 



•^2 1 / x(A M +, A p -) = ieu(qi, X^-)^v(q 2 , A M +) . (14) 

The part of the amplitude that comes from the final 
state emission can be written as 



M FSR (X e + , A e - , Ai) 



v\(pi,X e +)D u I (p 2 ,X e -) + vl I (pi,X e +)D + un(p2,X e -) 

(15) 

where the four-vector reads 



D»(\x) =ieF 2 „(s){ (qi + fa-q 2 ) 
+ (q2 + fa - q\) 



n gi ■ £*(fci, Ai) 



ijU q 2 ■ £*(fci, Ai) 
92 • fa 



qi ■ fa 
2e*'*(fci,Ai) } , (16) 



for 7r + 7r in the final state, while for in the final 

state it is 



D 11 (Ai , A„+ , A„- ) = ie 



u>i(q2, A A1 -)i M w/(gi, A M +) + u\ I (q 2 ,X fl -)B> J 'v II (q 1 ,X tl +) 



(17) 



with 



A" = 



(s*(fa, XJ-kt + 2e*(fci, Ai) ■ g 2 ) a*- 
2fa ■ q 2 

a"- (2e*{fa, Ai) ■ gi + k+e*{fa, AQ-) 
2fa ■ q 1 



(18) 



(e*(fci, Ai)+fcr + 2e*(fci, Ai) • g 2 ) ^ + 
2fci • g 2 

<t"+ (2e*(fci, Ai) ■ gi + k^e*{fa,X 1 )+) 
2fa ■ qi 



(19) 



The FSR matrix element squared and the FSR-ISR 
interference for pions in the final state agrees numerically 
(15 digits) with the code of EVA [7], if non-leading mass 
terms ~ m 2 missing in EVA are added. The largest rel- 
ative change of the matrix element squared due to those 
missing terms is however as small as 10~ 6 . The sum over 
polarisations of the squared matrix element for muon fi- 
nal states is numerically identical to the result obtained 
by means of the trace method using FORM 3|] . For both 
final states the external gauge invariance was checked nu- 
merically, while one can see at a glance that the above 
analytical formulae have that property. 

To analyse the contribution from FSR the program 
can be run in three different options: initial state radia- 
tion only, initial state radiation plus final state radiation 
without interference and complete result with interference 
terms. In the last two cases it is necessary to change the 
generation of the phase space to absorb final state emis- 
sion peaks. For these two options we use three channels 
to absorb the peaks in Q 2 and in pion (muon) angular 
distributions. For the muon case we use the approximant 



A(g 2 ,cos(^)) 

1 



1 l 1 

1 — q 2 q 2 
1 



1 - q 2 V 1 - "K) cos(6» M ) 1 + v(m ll ) cos(0 M ) 



where 



v(m) 



Am 2 



2_Q 2 



(20) 
(21) 



For the pion case we use 



f 2 {q 2 , cos(e n )) = 



1 



1 



2 (q 2 — m 2 / s) 2 + r 2 m 2 / s 2 
1 



1 — v(m 7r ) cos(#7r) 1 + w(m w ) cos^) 

(22) 

An appropriate change of variables allows for a smooth- 
ing of the aforementioned peaks. The q 2 and cos(0p( 7r )) are 
generated according to the functions /i(2)('Z 2 ,cos(0 Al ( Tr ))). 
These very simple approximants work well enough to al- 
low for relatively fast event generation. 



B The analytical formulae used in Section 2 

The formulae resulting from the analytical evaluation of 
the integration over photon angles are adopted from Refs. 
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29 3,0] (see also Section 2 for details). The contributions 
of the virtual + soft corrections to the hadronic invariant 
mass (Q 2 ) differential distribution are given by 



X ( 1 + - 

7T 



log(4 W 2 )(L-l) + ^-2 + ^ 



a 

17 



1 + q 4 



log(q 2 ) L 2 



1 + q 4 
l-q 2 



2(1 -q 2 ) 
(u 2 {l-q 2 ) + \og{q 2 )\og(l-q 2 ) 



\0g 2 (q 2 ) , 5 



+ -log(<? 2 ) -(l-^)log(^) + 



1 + q 4 



l-q 2 



S h2 (l-q 2 ) + 



log(g 2 



U 2 (l-q 2 ) 



log(aiog(l-? 2 )-^5^+M? 2 ) 



-31og(l-g ) - 8 — — 

+ (l + (? 2 )^2Li 3 (l-q 2 )-5i, 2 (l- ( Z 2 ) 

logV) 



log(l-^)Li 2 (l- (Z 2 ) + 



l-7g 2 
2 

1 - 5q 2 
4 

3-2g 2 



Li 2 (l-g 2 )+log( (? 2 )log(l- (Z 2 ) 

log 2 (l-<? 2 ) + ^ 
log(l-g 2 ) + I^!log( g 2 )-l 



(23) 



While the emission of two hard photons, i.e. both photons 
with energy larger than W\fs, contributes as 



{1 + q 2 ) (}ML r ^l + s lA i- q 2 ) 



TT 2 q 2 



U 2 (l-q 2 )-(10-25q^° g(q2) 



1 2q 2 

2 + 1T 

2 1 7q 2 

{l-q 2 ) 2 ~ 4 ~ T~ 
,2 o „2 



() 



logV) 



1-g 2 2 q 2 ( log(g 2 ) \ 2 - 
2 31-(7 2 ^ l- 9 2 ;. 



(24) 



with L = \og(s/m 2 ) and 51,2 the Nielsen's generalised 
polylogarithm function 

^W = t^/ I log»-'( t )^^*,(25) 
(n - 1)! pi J t 

Li„(z) = S n -i,i(z) being the polylogarithms 



Li, 



k=l 



(26) 



and 



S h2 {l-z) = -log 2 (z)log(l-z) + C(3) 
+ log(z)Li 2 (z) - Li 3 (z) . 



(27) 



The function R(Q 2 ) is related to the hadronic current 
jem through 



JfT(Jr)* d§ n (Q;q u ...,q n ) 



1 



^ (Q^Q* - 9^Q ) RiQ ) 



(28) 



where d<P n (Q; qi, . . . , q n ) denotes the n-body phase space 
with all statistical factors coming from the hadronic final 
state included. 

The ratio R(Q 2 ) = a(e + e~ — > hadrons)/cr po i„t for 
hadrons = 7r + 7r~ is equal to 



da 4a 3 
1 + q 4 



2lo g ( 1 -^)- 1 ^l\(L -I) 2 
1 — \ 8V 2w 2 r 

-H-(i- 9 2 ) + (i + ^^i^ 



+ < 1(1 - q 2 ) - q 2 \og{q 2 ) + (1 + 9 2 )^!l 
2 ~ ^1,2(1 - 9 ) 5 — Li 2 (l - q ) 



31og 2 ((? 2 ) ,7T 2 5v . , 2 , 



i?(Q 2 ) = |F(Q 



2\|2 



(29) 
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